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While graphene is a semi-metal, a recently synthesized hydrogenated graphene called graphane, 
is an insulator. We have probed the transformation of graphene upon hydrogenation to graphane 
within the framework of density functional theory. By analyzing the electronic structure for eighteen 
different hydrogen concentrations, we bring out some novel features of this transition. Our results 
show that the hydrogenation favours clustered configurations leading to the formation of compact 
islands. The analysis of the charge density and electron localization function (ELF) indicates that as 
hydrogen coverage increases the semi-metal turns into a metal showing a delocalized charge density, 
then transforms into an insulator. The metallic phase is spatially inhomogeneous in the sense, it 
contains the islands of insulating regions formed by hydrogenated carbon atoms and the metallic 
channels formed by contiguous bare carbon atoms. It turns out that it is possible to pattern the 
graphene sheet to tune the electronic structure. For example removal of hydrogen atoms along the 
diagonal of the unit cell yielding an armchair pattern at the edge gives rise to a band gap of 1.4 eV. 
We also show that a weak ferromagnetic state exists even for a large hydrogen coverage whenever 
there is a sublattice imbalance in presence of odd number of hydrogen atoms. 

PACS numbers: 73.22.Pr 61.48.Gh 81.05.ue 



I. INTRODUCTION 



Carbon is regarded as one of the most versatile ele- 
ments in the periodic table forming a wide variety of 
structures such as three dimensional sp 3 bonded solids 
like diamond, sp 2 hybridized two dimensional systems 
like graphene and novel nano structures like fullerenes 
and nanotubes. The electronic structure and the physi- 
cal properties of these carbon based materials are turning 
out to be exotic [l| . Although the existence and the prop- 
erties of three dimensional allotrope, graphite containing 
weakly coupled stacks of graphene layers were well known 
Q, the experimental realization of a monolayer graphene, 
brought forth a completely different set of novel proper- 
ties [3|, |4[ . The triangular bipartite lattice of graphene 
leads to the electronic structure having a linear disper- 
sion near the Dirac points. The low energy behaviour of 
such two dimensional electrons in graphene has been a 
subject of intense experimental and theoretical activity 
exploring the electronic, magnetic, mechanical, transport 
properties etc. For a recent review the reader is referred 
to Castro Neto et. al. [5|. 

Although graphene is considered as a prime candidate 
for many applications, the absence of a band gap is a 
worrisome feature for the applications to the solid state 
electronic devices. In the past, several routes have been 
proposed to open a band gap [6-8] . The most interesting 
one is the recent discovery of a completely hydrogenated 
graphene sheet named as graphane. Graphane was first 
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predicted by Sofo et. al. [9] on the basis of electronic 
structure calculations and has been recently synthesized 
by Elias et. al. [lOj . The experimental work also showed 
that the process of hydrogenation is reversible, making 
graphane a potential candidate for hydrogen storage sys- 
tems. Since upon hydrogenation, graphene, a semi-metal 
turns into an insulator, it is a good candidate for inves- 
tigating the nature of metal insulator transition (MIT). 

There are few reports investigating the electronic 
structure of graphene sheets as a function of hydrogen 
coverage leading to the opening of a bandgap. Most 
of these are restricted to a small number of hydrogen 
atoms. The electronic structure of hydrogen adsorbed 
on graphene has been investigated using density func- 
tional theory (DFT) by Boukhvalov et. al. [11] and Ca- 
solo et. al. [l2|. Their results support the use of 
graphene as a hydrogen storage material. Their work 
also shows that the thermodynamically and kinetically 
favoured structures are those that minimize the sublat- 
tice imbalance. Flores et. al. have investigated the 
role of hydrogen-frustration in graphane like structures 
using ab-initio methods and reactive classical molecular 
dynamics [l3j . The stability trends in small clusters of 
hydrogen on graphene have been discussed in a recent 
review [l4[. A recent work by Zhou et. al. [15] pre- 
dicted a new ordered ferromagnetic state obtained by 
removing hydrogen atoms from one side of the plane of 
graphane. Recently, Wu et. al. [l6| have reported impli- 
cations of selective hydrogenation by designing an array 
of triangular carbon domains separated by hydrogenated 
strips. The electronic structure of the interface between 
graphene and the hydrogenated part has been investi- 
gated by Schmidt and Loss [l7j. They show the exis- 
tence of edge states for zigzag interface having strong 
spin orbit interaction. An interesting experimental work 
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using; scanning tunnelling microscopy, by Hornekaer et. al 
fl8l . Tl9| demonstrates clustering of hydrogen atoms on 
graphite surface. The work shows that there is vanish- 
ingly small adsorption barrier for hydrogen in the vicinity 
of already adsorbed hydrogen atoms. The question of the 
effect of defect on the electronic structure of atomic hy- 
drogen has been addressed by Duplock et.al [20| within 
the density functional theory. They show that the elec- 
tronic gap state associated with the adsorbed hydro- 
gen is very sensitive to the presence of defect such as 
Stone- Wales defect. The localization behaviour of dis- 
ordered graphene by hydrogenation has been reported 
within the tight binding formalism [2l[ . The potential of 
hydrogenated graphene nanoribbons for spintronics ap- 
plications has been investigated within DFT, by Soriano 
et. al [22j . The defect and disorder induced magnetism in 
graphene (including adsorbed hydrogen as a defect) has 
been studied by Yazyev, Yazyev and Helm mostly using 
Hubbard Hamiltonian and DFT [23l-l25|. Theoretical in- 
vestigations have also been carried out for a single hydro- 
gen defect on graphane sheet using GW method [26] , and 
for one and two vacancies in graphane [27j using DFT. 
It has been shown that interaction between adatoms in 
hydrogenated graphene is long range and its nature is 
dependant on which sublattice the adatoms reside (28| . 
The result suggests that the adatoms tend to aggregate. 
A very recent work explores a formation of quantum dots 
as small island of graphene in graphane host [29^ 

In the present work, we investigate some aspects of 
graphene - graphane transition by probing the electronic 
structure of hydrogenated graphene. The objective of 
our work is to understand the evolution of the electronic 
structure upon hydrogenation and gain some insight into 
the way band gap opens. Therefore we have carried 
out extensive calculations for eighteen different hydrogen 
coverages between graphene (0% hydrogen coverage) and 
graphane (100% hydrogen coverage) within the frame- 
work of DFT. Our results suggest that the hydrogena- 
tion in graphene takes place via clustering of hydrogens. 
Analysis based on density of states (DOS) indicates that 
before the gap opens (as a function of hydrogen cover- 
age) the hydrogenated graphene sheet acquires a metal- 
lic character. This metallic state is spatially inhomoge- 
neous in the sense, it consists of insulating regions of 
hydrogenated carbon atoms surrounded by the conduct- 
ing channels formed by the bare carbon atoms. We also 
show that it is possible to tune the electronic structure 
by the selective decoration of hydrogen atoms to achieve 
semiconducting or metallic state. 

In the next section we present the computational de- 
tails, followed by results and discussions. 



II. COMPUTATIONAL DETAILS 

The calculations have been performed using a plane- 
wave proj ector augmented wave method based code, 
VASP [30[. The generalized gradient approximation as 



TABLE I: Shifted binding energies (BE) for different config- 
uration of 20% coverage (10H + 50C). The binding energy 
of the most stable structure is the reference (zero) level. The 
more negative binding energy means more unstable structure. 



20% Hydrogen Coverage 


Configuration 


BE (total cell) (eV) 


Random placement of hydrogens 


-9.32 


Five hydrogen pairs placed 
separately 


-4.29 


Hydrogens placed along the diago- 
nal of the cell 


-2.85 


Two clusters of seven and three 
hydrogens 


-2.15 


A single hydrogen isolated from the 
cluster of nine hydrogens 


-1.77 


One compact cluster 


0.0 



proposed by Perdew, Burke and Ernzerhof [31] has been 
used for the exchange-correlation potential. The con- 
vergence of binding energies with respect to the size of 
the super cell has been checked using three different sizes, 
viz., 5 x 5, 6 x 6 and 7x7 containing 50, 72 and 98 car- 
bon atoms for the case of 50% hydrogen coverage. The 
binding energy per atom changes by about 0.05 eV (a 
percentage change of 0.06%) in going from 5 x 5 to 7 
x 7. Therefore we have chosen a 5 x 5 unit cell for the 
coverage upto 50% of hydrogen and 6 x 6 for the higher 
coverages. This choice is consistent with the one used by 
Lebegue et.al [26]. In order to obtain adequate conver- 
gence in the density of states, we have carried out cal- 
culations on different k grids. It was found that at least 
9 x 9 k grid was required during geometry optimization 
for an acceptable convergence. However, a minimum of 
17 x 17 k grid was necessary for obtaining an accurate 
DOS. The convergence criterion used for the total energy 
and the force are 10 -5 eV and 0.005 eV/A respectively. 

All the calculations have been performed on the chair 
conformer configuration where hydrogen atoms are at- 
tached to carbon atoms alternatively on opposite sides 
of the plane. This is known to be a lower energy config- 
uration as compared to the boat conformer [9|, [26| . 

III. RESULTS AND DISCUSSION 

In order to decide the minimum energy positions for 
hydrogen atoms, following procedure has been adopted. 
Upto 20% coverage of hydrogen, we have carried out 
the geometry optimization for two different configura- 
tions of hydrogen atoms, first by placing the hydrogen 
atoms randomly and second by placing the hydrogen 
atoms contiguously, so as to form a compact cluster of 
hydrogenated carbon atoms. It turns out that, in all the 
cases, the configuration forming the compact cluster of 
the hydrogen atoms is energetically favoured. In order to 
assess the relative stability of different patterns we have 
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TABLE II: Shifted binding energies (BE) for different con- 
figuration of 50% coverage as shown in figure U The binding 
energy of the most stable structure is the reference (zero) 
level. The more negative binding energy means more unsta- 
ble structure. 



50% Hydrogen Coverage on 98 Carbon cell 


Configuration 


BE (total cell) (eV) 


Random placement of hydrogens 
(figure [Ha)) 


-23.29 


Cham or connecting bare Carbon 
atoms (figure E^b)) 


-3.65 


Three separated islands (figure 

Wi) 


-3.25 


Zigzag interface with a line of bare 
carbon atoms at the two side edges 
(figure TLd)) 


-2.46 


Compact mixed (zigzag & arm- 
chair) interface (figure H{e)) 


-1.66 


Compact armchair interface (figure 

mo) 


-1.43 


Compact zigzag interface (figure 

Kg)) 


0.0 



calculated six different patterns for 20% coverage case as 
shown in table HI The binding energy of the most sta- 
ble structure is the reference level. The more negative 
binding energy means more unstable structure. 

We note that the compact cluster configuration is the 
lowest in energy and is lower by 0.9 eV/H atom compared 
to the energy of the configuration with randomly placed 
hydrogen atoms. Thus our geometry optimization shows 
a preference for hydrogens to decorate the graphene lat- 
tice in a contiguous and compact manner. In order to 
assess the validity of this process for the larger coverages 
and to understand the influence of different edge patterns 
like zigzag and arm chair, we have considered seven dif- 
ferent patterns for 50% coverage, as shown in figure [TJ 
These calculations have been carried out on a larger unit 
cell containing 98 atoms of carbon. Although a compact 
cluster configuration is energetically preferred one, there 
are different ways in which the hydrogens can be placed 
yielding a compact geometry. The seven cases considered 
are : 1) Randomly distributed hydrogens (figure a)), 
2) A chain of bare carbon atoms (figure QJb)), 3) Three 
separated clusters as shown in (figure QJc)), 4) A single 
cluster having zigzag interface and with one line of bare 
carbon atoms at the two side edges (figure DJd)), 5) A 
mixture of armchair and zigzag pattern at the interface 
(figure Die)), 6) A armchair pattern at the interface (fig- 
ure [Hf)) and 7) A zigzag pattern at the interface (figure 

In the zigzag case each hexagonal ring at the edge has 
three hydrogenated carbons and three empty carbons. 
An armchair pattern consists of only two hydrogenated 
carbons at the edge. The binding energies for all these 
seven structures (shown in figure [1]) are compared in the 



TABLE III: Shifted binding energies (BE) (eV) and the 
magnetic moments(/x j g) of 4 hydrogen vacancies. The binding 
energy of the most stable structure is the reference(zero) 
level. The more negative binding energy means more unsta- 
ble structure. 



92% Hydrogen Coverage (4 vacancies) 


Configuration 


BE (total cell) 


Magnetic Moment 


Random 


-5.50 


4.0 


3H on one hexagon 


-2.09 


1.82 


and 1 on other 






2H pairs placed on 
different hexagons 


-0.11 


nil 


Compact single 
island 


0.0 


nil 



table HIl Even among the islands, the 'most compact' 
one (having minimum surface area of covered hydrogens) 
has the highest binding energy. By the most compact, we 
mean an island of hydrogenated carbons which covers the 
minimum area. We have discussed more details about the 
zigzag and arm chair patterned configurations in section 

Em 

The tendency to form the compact cluster can be un- 
derstood by noting that the hydrogen atom placed on 
the top of an bare carbon atom pulls up the carbon atom 
above the plane by 0.33 A and deforms the surrounding 
lattice points also. Therefore it costs less energy to place 
an extra hydrogen atom nearest to the existing cluster, 
since the number of strained bonds is less compared to 
the case where the hydrogen atom is placed away from 
the cluster. In the latter case the entire neighbourhood 
is deformed. Our results are consistent with the work of 
Hornekaer et al [18] . 

We have also carried out a similar analysis for the case 
of four hydrogen vacancies in graphane. The calculations 
for various configurations have been carried out with spin 
polarization. The calculated binding energies and the 
magnetic moments are tabulated in table IIII1 The four 
configurations considered are 1) 4 hydrogen atoms re- 
moved randomly, 2) 3 hydrogen atoms removed from one 
hexagon and 1 from other, 3) 2 hydrogen atom pairs re- 
moved from different hexagons and 4) 4 hydrogens re- 
moved from a single hexagon (compact). We observe 
that the structure having a compact form of vacancies 
has the highest binding energy. The magnetic moment 
is about 2 for the lattice imbalance case. It is fruitful 
to recall Lieb's theorem which states that for a bipartite 
lattice (one electron per site) the spin S of the ground 
state is 1/2 x (the lattice imbalance) [32j. The mag- 
netic moment seen in the case of random placement is 
just the sum of the isolated moments of hydrogen atoms, 
essentially a non interacting case. 

Now we discuss the total DOS for different hydrogen 
coverages which are shown in figure [2] and [3j All the 
DOS are for the non spin polarised cases. The geometry 
used is for the minimum energy configuration (compact 
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FIG. 1: (Colour online) Different hydrogen decorations for 50% concentration. In figure, yellow (in print light shaded) balls 
are bare carbon atoms, turquoise (in print darker shades) balls are hydrogenated carbons and red (in print dark small) balls 
are hydrogen atoms. 



cluster of hydrogen atoms). In both the figures the plots 
are shown in a restricted region to enhance the clarity 
near the Fermi energy. The effect of addition of a small 
concentration of hydrogen in the V-shape DOS near 
can be seen from figure (0(b) toEfd)). It can be noted 
that the characteristic V shape valley seen in figure [2fa) 
is due to peculiar sp 2 bonded carbon atoms in graphene, 
and the addition of hydrogen atoms immediately distorts 
the symmetry producing the deformation in the DOS. 
Thus in this region, the DOS is modified by additional 
localized a-p z bonds between the hydrogen and carbon. 

As the hydrogen coverage increases there is a signifi- 



cant increase in the value of DOS at the Fermi level. The 
process of hydrogenation is accompanied by the change in 
the geometry. The hydrogenated carbon atoms are now 
moved out of the graphene plane, in turn the lattice is 
distorted and the symmetry is broken. As a consequence, 
more and more k points in the Brillouin zone contribute 
to the DOS near Fermi level, the increase being rather 
sharp after 20% coverage. The region ranging from 30% 
coverage to about 60% coverage is characterized by the 
finite DOS of the order of 2.5/eV near the Fermi energy. 
As we shall discuss this region can be described as having 
metallic character with delocalized charge density. 
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FIG. 2: (Colour online) The total DOS for hydrogenated graphene for various hydrogen concentrations below 90%. The zero 
of the energy is taken at the Fermi level and is marked by a vertical line. X axis denotes E — Ef. 
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FIG. 3: (Colour online) The total DOS for hydrogenated graphene for various hydrogen concentrations above 90%. The zero 
of the energy is taken at the Fermi level and is marked by a vertical line. X axis denotes E — Ef. 



Now we complete the discussion of DOS by presenting 
the cases of hydrogen coverages greater than 90%. In 
figure [3] we show the density of states obtained by the re- 
moval of one - four hydrogen atoms in a unit cell contain- 
ing fifty carbon atoms. Clearly one and three hydrogen 
vacancies (i.e 98% hydrogen and 94% hydrogen respec- 
tively) induce states on the Fermi level while two and four 
hydrogen vacancies (i.e 96% hydrogen and 92% hydrogen 
respectively) do not induce any state at the Fermi level 
(See the discussion of hydrogen imbalance later). For 
a small number of vacancies these are 7r bonded states 
localized on bare carbon atoms. 

In summary, it is possible to discern, rather broadly, 
three regions of hydrogen coverage. A low concentration 
region, where the DOS undergo the distortions due to 
loss of lattice symmetry of the system, intermediate con- 
centration region where a metallic-like phase is seen and 
a vary high concentration region where most of the car- 
bon atoms are hydrogenated and vacancy gives rise to 
midgap states. 

As discussed before, our calculations show the pres- 
ence of localized states in the DOS for a very low hydro- 



gen coverage as well as for very high hydrogen coverage 
(midgap states). These DOS show a pattern of peak or 
a valley at Ef (in case of low coverage) and at the centre 
of the gap (in case of high coverage). The reasons can 
be attributed to the existence of sublattice imbalance. If 
the difference of hydrogen atoms on the two sides of the 
sheet is odd then it leads to a peak. Hence by adding 
impurities one by one to graphane we get a sequence of 
midgap states. This is in consistent with the work of Ca- 
solo et al [12| . Interestingly this feature is also retained 
for the intermediate coverages of hydrogen (where we get 
a finite DOS at the Fermi level) when the sub lattice im- 
balance of hydrogen is one, e.g., in figure Efe), [2ff) and 
Etg)? we see the peaks at Fermi level. As we shall see 
in section IIII Al this feature has an implication for the 
magnetic behaviour of the system. 

In order to bring out the difference between the local 
electronic structure of the bare carbon atoms and the 
hydrogenated carbon atoms, we have analyzed the site 
projected DOS for all the cases. In figure HJ we show 
site projected DOS for 40% hydrogen coverage depict- 
ing the contributions from a hydrogenated carbon site 
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FIG. 4: (Colour online) Site projected DOS for hydrogenated 
carbon sites (dotted line) and bare carbon sites (continuous 
line) for hydrogen coverage of 40%. Almost all the contribu- 
tion comes from the p z which has been showed in the figure. 
Note that only bare carbon atoms contribute to the DOS at 
Fermi level. 




20 40 60 80 100 
% of Hydrogen 

FIG. 5: (Colour online) Variation of the value of DOS at the 
Fermi level as a function of hydrogen coverage. 



and a bare carbon site. Quite clearly, the contribution 
around the Fermi level comes from the p z orbitals of bare 
carbon atoms only. This is a general feature for all the 
systems investigated. It may be emphasized that in pure 
graphene all the carbon atoms contribute to a single k 
point (Dirac point). In contrast to this, upon hydrogena- 
tion, only bare carbon atoms contribute to the DOS at 
Fermi energy, and they do so at many k points of the 
Brillouin zone. As the concentration increases further 
(above 80%), there are too few bare carbon atoms avail- 
able for the formation of delocalized 7r bonds. The value 
of DOS approaches zero and a gap is established with a 
few mid gap states. 

The evolution towards the metallic state can be better 
appreciated by examining the variation in the value of 
DOS at the Fermi level which is shown in figure [5] A clear 
rise in DOS is seen after 20% hydrogen concentration 
peaking around 3.5 eV at 50% concentration. This rise 
is due to the increasing number of k points contributing 



to the Fermi level, as inferred from the analysis of the 
individual bands. Evidently, over a significant range of 
concentration, the value at Fermi level is more than 2/eV. 
The decline seen after 60% is because of the reduction in 
the number of bare carbon atoms. The character of the 
DOS changes after about 80%. The value of the DOS 
at the Fermi level oscillates between zero and some finite 
value. It is most convenient to describe this region as 
graphane with defects by the removal of a few hydrogen 
atoms giving rise to mid gap states. The nature and 
the placement of the induced states is dependent on the 
number of hydrogen atoms removed. 

It may be emphasized that the presence of states 
around the Fermi level giving finite DOS does not guar- 
antee that the system is metallic unless we examine the 
nature of localization of the individual states. Therefore 
we have examined the energy resolved charge densities of 
the states near the Fermi level and the electron localiza- 
tion function (ELF). The energy resolved charge densities 
are obtained by summing the charge densities of all the 
k points contributing in the small energy region near Ef. 
Therefore, the charge densities shown in (figure [6j) for 
40%, 70% and 92% hydrogen coverages specifically brings 
out the nature of the states near Ef. A particularly strik- 
ing feature is the formation of two spatially separated re- 
gions as seen in figure [6^b) and[6^c). The hydrogenated 
regions hardly contribute to the charge density giving 
rise to the insulating regions surrounded by the 7r bonded 
bare carbon atoms forming conducting regions. This fea- 
ture is also seen for the higher concentrations upto 70%. 
It may be emphasized that the topology in this range 
of concentrations (30%-70%) shares a common feature 
namely, there is a contiguous region formed by the bare 
carbon atoms. It may be pointed out that the contigu- 
ous charge density is attributed to the favoured configu- 
ration of compact cluster formation [38|. The change in 
the character of the state at 90% and above is also evi- 
dent in figure [6jd). There are insufficient number of bare 
carbon atoms to form contiguous regions. As a conse- 
quence, these carbon atoms form localized bonds giving 
rise to mid gap states noted earlier. 

The degree of derealization of an electron of the sys- 
tem can be understood by examining the electron local- 
ization function as follows. For a single determinant al 
wave function built from Kohn-Sham orbitals, the 
ELF is defined as, [H 

XELF = [l + (^/^) 2 ]- 1 , (1) 

where 

D h = (3/10)(3tt 2 ) 5/ V /3 , (2) 

D = (1/2) £ |V^| 2 - (1/8)| Vp\ 2 /p, (3) 

with p = p(r) being the electron density. D is the 
excess local kinetic energy density due to Pauli repulsion 
and D/,, is the Thomas-Fermi kinetic energy density. The 
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FIG. 6: (Colour online) Isosurfaces of charge densities of bands near the Fermi level (see text). For comparison, the charge 
density of graphene at Dirac point is also shown. 




FIG. 7: (Colour online) ELF plots with two different isosurface values for 70% hydrogen concentration. Figure (a) shows for 
isosurface value 0.75 and figure (b) shows for isosurface value 0.62 . 



numerical values of Xelf are conveniently normalized to 
a value between zero and unity. A value of 1 represents a 
perfect localization of the charge, while the value for the 
uniform electron gas is 0.5. Typically, the existence of an 
isosurface of a high value of Xelf say, > 0.70, signifies a 
localized bond. 

We have examined the isosurfaces of ELF for various 
values between 0.5 and 1.0. The isosurface for a typi- 
cal value 0.75, is shown in figure [7fa) for 70% coverage. 
It can be seen that for this high value isosurface exists 
mainly along carbon- hydrogen bonds (a-p z ). The exis- 
tence of localized a bonds between the bare carbon atoms 
can also be noticed. In figure Efb) we show the isosur- 
face for the value of 0.62. The figure shows a continuous 



surface covering all the bare carbon atoms. This signifies 
the delocalized nature of the charge density arising out 
of 7r bonds formed by p z orbit als. We have analyzed the 
ELF for all the coverages. The above features are seen 
for the coverages from 30% to 80%. Thus the analysis of 
ELF confirms the delocalized nature of the charge den- 
sity near the Fermi level for the concentration from 30% 
to 80%. 



A. Magnetism 

The existence of a peak at Fermi level in the non spin 
polarized calculation is an indication of Stoner instability 
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FIG. 8: (Colour online) (a) Spin density, (b) total DOS, (c) and (d) site projected (p z only) DOS for two bare carbon atoms 
from two different sublattices. X axis of the plots denotes E — Ef. This figure is for 50% hydrogen coverage. 



that may lead to a more stable spin polarized solution. 
As an example we have analyzed the case of 50% hy- 
drogen coverage which shows a peak in the DOS at the 
Fermi level. When we allow spin polarization in our cal- 
culations, indeed we find a stable ferromagnetic solution, 
e.g., for 50% hydrogen coverage we get a total magnetic 
moment ~ 1 /ie per unit cell as reflected in the spin den- 
sity plot shown in figure [5Ja). However, the exchange 
splitting seen on the bare carbon atoms is rather small 
and is expected to survive only at a very low temperature. 
This fragile nature of magnetism is clearly indicated by 
the collapse of magnetic moment with a smearing width 
of 0.08 eV. The site projected DOS for two bare carbon 
atoms belonging to two different sublattices along with 
the spin density plots are shown in figure (|8jc) and[8^d)). 
Carbon atom belonging to one sublattice shows a spin- 
polarized behaviour whereas the other sublattice carbon 
atom has an energy gap in the electronic spectrum. This 
sublattice effect is similar to what is observed in case of 
graphene nanoribbons [34j . 



B. Tuning the electronic structure with 
hydrogenation 



Now we bring in another interesting aspect of hydro- 
genation brought out by our calculations. Our results 
indicate that it is possible to pattern the graphene lat- 
tice with hydrogen and tune the electronic structure. In 
figure [9ja), the patterning is done by removing the hy- 
drogen atoms along the diagonal of the unit cell where 
as in figure [9f b) the hydrogen atoms are removed paral- 
lel to one edge of the unit cell. The horizontal pattern 
is more stable than the diagonal pattern by 0.04 eV/H 
atom. The corresponding DOS are shown in figure [9fc) 
where a dramatic difference is seen. The diagonal pat- 
tern shows a clear band gap of 1.4 eV whereas the hor- 
izontal pattern shows finite DOS along with a magnetic 
solution. A close inspection of the geometry reveals an 
analogy with graphene nano ribbons (GNR). The diag- 
onal pattern resembles an armchair GNR with a chain 
of hexagons while the horizontal pattern is analogous 
to a zigzag GNR with a width of 3 rows in this par- 
ticular case. The corresponding DOS also resemble the 
electronic structure of an armchair (zigzag) GNR with a 
semiconducting (metallic) nature j35j. It is a reasonable 
conjecture that the band gaps of the patterned system 
can be tuned by controlling the width of the bare carbon 
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FIG. 9: (Colour online) Decoration of hydrogen (a) along the diagonal of the unit cell and (b) along the edge of the unit cell. 
In figure, yellow (in print light shaded) balls are bare carbon atoms, turquoise (in print darker shades) balls are hydrogenated 
carbons and red (in print dark small) balls are hydrogen atoms, (c) DOS plots for the two cases. The lone atoms seen are from 
repeating super cells 



channels analogous to the case of armchair GNR where 
the band gap decreases as the width of GNR is increased. 



IV. CONCLUDING REMARKS 

In conclusion, our detailed density functional investi- 
gations have revealed some novel features of graphene- 
graphane metal-insulator transition. As the hydrogen 
coverage increases, graphene with a semi metallic charac- 
ter turns first into a metal and then to an insulator. Hy- 
drogenation of graphene pulls the carbon atom out of the 
plane breaking the symmetry of pure graphene. As a con- 
sequence, many k points in the Brillouin zone contribute 
to the DOS at the Fermi level giving rise to a metallic 
system. The metallic phase has some unusual character- 
istics: the sheet shows two distinct regions, a conducting 
region formed by bare carbon atoms and embedded into 
this region are the non conducting islands formed by the 
hydrogenated carbon atoms. The onslaught of insulating 
state occurs when there are insufficient numbers of bare 
carbon atoms to form connecting channels. This also 



means that the transition to insulating phase depends on 
the distribution of hydrogen atoms and will occur when 
the continuous channels are absent. The present work 
opens up the possibility of using partially hydrogenated 
graphene having designed patterns of conducting chan- 
nels along with insulating barriers for the purpose of de- 
vices. Our results also show that it is possible to design a 
pattern of hydrogenation so as to yield a semiconducting 
sheet with a band gap much lower than that of grapheme. 
Finally we may note that the calculation of conductiv- 
ity in such a disordered system is a complex issue. The 
present study focusses on the evolution of the density of 
states to understand the change in the character of single 
particle orbit als as a function of hydrogen coverage. An 
obvious extension of this work is the study of transport 
properties to have a more vivid picture. 
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